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Abstract. We present the estimate of the evolution of type 1 AGNs in 
the hard (2-10 keV) X-rays drawn from the HELLAS survey, and in the IR 
(15 //m) obtained from the ELAIS survey. We find that the local luminos- 
ity function (LF) of AGNl in the 2-10 keV is fairly well represented by a 
double-power-law-function. There is evidence for significant cosmological 
evolution according to a pure luminosity evolution model -Lx(-z)oc(l + 
z)'', with A:=2.12t°;u and k=2.19toil in a (l^m,f^A)=(1.0,0.0) and in 
a (nm5^A)=(0.3,0.7) cosmology respectively. In a (r2nnf^A)=(l-0,0.0) 
Universe the data show an excess of faint high redshift type 1 AGN 
which is well modeled by a luminosity dependent density evolution, sim- 
ilarly to what observed in the soft X-rays. In the IR band, with a 
(f2ni)f^A)=(l-0,0.0) cosmology, the evolution is similar to what observed 
at other wavebands, the LF is a double-power-law-function with a bright 
slope 2.9 and a faint slope 1.1, following a pure luminosity evolution model 

+0.16 
-0.20- 



with A:=3.00+°-^^ 



1. Introduction 

Before the advent of the last generation of hard X-ray telescopes, AGN samples 
were predominantly based on type 1 objects (AGNlQ) selected either in the 



^HELLAS: A. Antonelli, A. Comastri, P. Giommi, G. Matt, S. Molendi, G.C. Perola, F. Pom- 
pilio. ELAIS: C. Lari, M. Mignoli, G. Zamorani, D. Alexander, S. Oliver, S. Serjeant, A. 
Franceschini, G. Danese, M. Rowan-Robinson. 

^We here define AGNl all the AGNs having a broad emission line spectrum in the optical. 
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optical or later on in the soft X-rays by ROSAT. In these bands the evolution 
of AGNl has been well measured (see e.g. Delia Ceca et al. 1992; Boyle et 
al. 2000; Miyaji, Hasinger and Schmidt 2000). On the contrary the generation 
of samples of type 2 AGNs (AGN2) has been difficult at any wavelength and 
limited to few local surveys. 

The general picture was in favor of a model in which type 1 objects were 
associated to AGNs with low absorption in the X-rays while type 2 to obscured 
sources with large column densities and spectra strongly depressed in the soft 
X-rays and emitting mainly in the a) hard X-rays and b) mid and far infrared. 
For these reasons, in order to achieve a more exhaustive picture of the evolution 
of the whole family of AGNs, the researches have been recently concentrated in 
these two wavebands. 

Thanks to their high angular resolution (~l-5 ") the first spectroscopic 
identification projects of faint Chandra and XMM-Newton hard X-ray sources 
have been able to observe faint (I~23) optical counterparts. At variance with 
the classical type-l/type-2 model in the optical, a significant number of the 
counterparts (~30%) resulted to be apparently optical normal galaxies, and 
moreover part of the optical type 1 AGNs resulted to be absorbed in the hard 
X-rays (see e.g. Fiore et al. 2000; Comastri et al. 2001; Barger et al. 2001). 

These observations have complicated the picture of the AGN model. In 
the framework of the computation of the density of AGNs, it is not clear how 
to handle the classification of the sources and to take into account the selection 
biases introduced by the observation in the 2-10 keV range where the absorption 
still play a relevant role. 

At variance, the situation in the mid infrared is in a less advanced stage 
because of the scarcity of satellite missions, the only ones being IRAS and ISO. 
Although the ISO mission has been completed a couple of years ago, the only 
existing survey of AGN in the mid infrared is still the sample from Rush, Malkan 
and Spinoglio (RMS, 1993) obtained from the IRAS observations of the whole 
sky at 12 /xm down to ~300 mJy. We have thus decided to undertake two 
projects. 

a) While the recent deep surveys with Chandra and XMM-Newton have 
reached fluxes ~ 2 x lQ~^^erg cm'^'^s^^ (2-8 keV, Brandt et al. 2001) in quite 
small areas (less than 1 deg^), at brighter fluxes (~ W~^^erg cm~^s~^; 5-10 
keV) the density of sources is low (about 5/deg^) and tenths of square degrees 
are to be covered in order to collect statistical significant samples. Such samples 
are quite rare at moment but are essential in order to fully cover the X-ray 
luminosity /redshift plane and study the AGN evolution. In the first part of 
this paper we report the results of the spectroscopic identifications of one of 
such samples (HELLAS), which covers from 0.2 to 55 deg^ in the flux range 
5 X 10-^^ - 5 X IQ-^^erg cm-^g-i (5-10 keV). 

b) In order to extend at deeper fluxes the coverage of the luminosity redshift 
plane in the mid infrared (15 /im) we have carried out the European Large Area 
ISO Survey (ELAIS) over 12 deg^ down to ~0.5-l mJy, more than 2 order of 
magnitudes deeper than the previous RMS sample. In the second part of this 
paper we will report on the results of the spectroscopic identifications of this 
project. 
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2. HELLAS: The evolution of AGNl in the Hard X-ray 

2.1. The HELLAS survey 

The X-ray sources have been detected by the BeppoSAX-MECS instruments in 
the 5-10 keV band in the framework of the High Energy Large Area Survey 
(HELLAS). Prehminary results have been presented in Fiore et al. (1999). The 
whole survey and the catalogue used in this study is described by Fiore et al. 
(2001). The data have been analyzed in the framework of the synthesis models 
for the X-ray background by Comastri et al. (2001), and the correlation with 
the soft X-rays has been investigated by Vignali et al. (2001). 

The spectroscopic follow up of the HELLAS sources has been carried out 
in a subsamplc of 118 sources out of a total of 147. The BeppoSAX X-ray 
positions have an uncertainty of about 1-1.5 arcmin depending on the distance 
of the sources from the axis of the telescope. We have thus searched for optical 
counterparts having R magnitude brighter than 20.5-21.0 in a circular region of 
1-1.5 arcmin of radius around the HELLAS positions. In the case of large off-axis 
distances, the larger error-boxes (1.5") have been used. 25 sources have been 
identified with cross-correlation with existing catalogues available from NED, 
and 49 have been investigated at the telescope. 

We have divided the identified sources in AGNl, AGN2, galaxies with nar- 
row emission lines showing modcratc-to-high degree of excitation (ELG), and 
stars. AGN2 includes AGN1.5, AGN1.8 and AGN1.9. No distinction has been 
done among radio-loud and radio quite objects. The R-band magnitudes of 
these identifications have been pushed down to the limit of obtaining less than 
10% probability of having a spurious identification due to chance coincidences 
for each of the HELLAS sources (i.e. 90% reliability for the whole sample). 

For AGNl wc have chosen a limit of R=21 where the surface density is 
~60-70/deg^. Wc have chosen a limit of R=19 for AGN2. We have collected 
all the objects which are mainly starburst galaxies and low ionization narrow 
emission line galaxies (LINERS) within the class of emission line galaxies (ELG). 
For these objects a limit of R=17.5 has been chosen. In total, 63% (74/118) of 
our HELLAS subsample has been searched for spectroscopic identification. 61 
have been identified: 37 AGNl, 9 AGN2, 5 ELG, 6 Clusters, 2 BL Lac, 1 Radio 
Galaxy and 1 Star. 

2.2. The evolution of AGNl 

For the reasons described before we have decided, at moment, to limit the anal- 
ysis of the evolution of AGN to the AGNl only. We have assumed Ho=50 
km/s/Mpc and either the (rim,^^A)=(l-0,0.0) or the (r2nij^^A)=(0.3,0.7) cosmolo- 
gies. In order to join our data with other existing samples we have computed 
the 2-10 keV luminosity for type 1 AGN assuming a typical spectral energy dis- 
tribution as computed by Pompilio, La Franca &; Matt (2000), which is roughly 
approximated by a single power slope in energy with index 0.6 (^ oc E~^'^) in 
the range 2-50 keV (i.e. up to 2=4). 

The statistical significance of our analysis has been increased by joining our 
data to other hard X-ray samples of AGNl identified by Grossan (1992), Boyle 
et al. (1998), and Akiyama et al. (2000). 
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We carried out a least-squares method to derive best-fit parameters for the 
2-10 keV luminosity function (LF) and its cosmological evolution. All the fits 
have been tested with the bidimensional Kolmogorov-Smirnov test (2DKS). The 
summary of the results is presented in Table 1. Two different functional forms 
have been used. Following Boyle et al. (1998) and Ceballos and Barcons (1996), 
we used the pure luminosity evolution (PLE) for the QSO LF. We also used the 
luminosity dependent density evolution (LDDE) model similar to the one fitted 
in the soft X-rays by Miyaji et al. (2000). 

The local {z=0) QSO LF used for the PLE model has been represented by 
a two-power-law: 

^ = A(Lj^-^^))L7^ ■.Lx<L*^{z = 0), 
'^ = AL-^' ■.Lx>L*^{z = 0), 

where Lx is expressed in units of lO^^ergs"^. A standard power-law luminosity 
evolution model was used to parameterize the cosmological evolution of this LF: 
L*^{z) = L*^mi + z)K 

In the case of the LDDE model, as an analytical expression of the present 
day {z = 0) luminosity function, we used the smoothly connected two power-law 
form: 

The description of the evolution is given by a factor e(Lx, z) such that: 
d^{L^,z) d^L^,z = 0) 



c/LogLx c/LogLx 



X e{L^,z), 



where 



e{L^,z) = 

(1 + ^)max(0,pl-aLog [La/Lx]) < < 

(1 + ^)^^ {z<z,;L^>L^) 
e(Lx,^c)[(l + ^)/(l + ^c)f' {z>z^). 

We started our computation of the luminosity function in an (r2m,r2A) = (l-0,0.0) 
cosmology. The only previous estimate of the evolution of AGNl in the 2-10 keV 
range was from Boyle et al. (1998), who joined the local sample of 84 AGNl 
observed by HEAOl from Grossan (1992) with a fainter sample of 12 AGNl 
observed by ASCA (see their distribution in the hx-z plane in Figure la). 

In Figure lb the luminosity function from only our 37 AGNl from HELLAS 
in three redshift intervals (0.0 < z < 0.2, 0.2 < z < 1.0, and 1.0 < 2; < 3.0) is 
shown. For comparison the best fit luminosity function from Boyle et al. (1998) 
is also shown. The data have been represented by correcting for evolution within 
the redshift bins as explained in La Franca and Cristiani (1997). The data are 
in rough agreement with the previous estimate of Boyle et al. (1998), but show 
some evidence of an excess of faint AGNl at redshifts larger than z ~0.3. This 
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Figure 1. a) The luminosity/redshift distribution for all the 156 
AGNl used to compute the hard X-ray 2-10 keV luminosity function. 
They include 37 AGNl from HELLAS, 82 AGNl from Grossan, 12 
AGNl from Boyle ct al. (1998) and 25 AGNl from Akiyama et al. 
(2000). b) The luminosity function of the 37 type 1 AGN from HEL- 
LAS. The dotted lines show, as a reference, the previously estimated 
luminosity function in the 2-10 kcV from Boyle et al. (1998). The 
densities have been corrected for evolution within the redshift bins. 



feature is strengthened by joining our data with the other AGNl samples from 
Grossan (1992), Boyle et al. (1998), and Akiyama et al. (2000) (see Figure 
2a). These samples all together collect 156 AGNl. These data have a (x^) 0.04 
probability to be drawn from a PLE model such as that computed by Boyle et al. 
(1998, see model in Table 1). Our best fit to the whole data with a PLE model 
found similar parameters as those of Boyle et al. (1998) but a slightly larger 
evolution {k = 2.13 instead of A; = 2.04) and a significantly larger normalization 
(model 1). This difference is originated by the necessity of better fitting the 
observed higher density of faint AGNl at high redshifts. The 2DKS test gives 
a probability of 0.22 for this fit (sec Figure 2a and Tabic 1). A even better 
probability of 0.31 is obtained if a stop in the evolution is applied at redshift 
Zcut = 1-39 and a larger evolution {k = 2.52) is used (model 2). 

Although our fits are already statistically adequate, our PLE model is not 
able to fully describe the observed over-density of faint AGNl at high redshifts. 
As this feature is quite similar with what observed in the soft X-rays (Miyaji 
et al. 2000) we tried to get a even better fit by using the luminosity dependent 
density evolution (LDDE) model similar to the one fitted in the soft X-rays. We 
first kept fixed all the parameters to the values found in the soft X-rays, just 
leaving free to vary (model 3). The 2DKS test gives a probability of 0.11. The 
value LogL*=44.13 found by our fit nearly corresponds to the value LogL*=43.78 
found in the soft X-rays if a hnear relation Lx(0.5-2 keV)=0.471xLx(2-10 keV) 
is assumed, which corresponds to the slope 0.6 of the AGNl spectra which we 
used in this analysis. Later we kept fixed just Zc = 1.55 and a = 2.5 to the 
value found by Miyaji et al. (2000) for AGNl only, and left all the remaining 
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Figure 2. a) The luminosity function of tlic total sample of 156 
AGNl including the 37 AGNl from HELLAS, 82 AGNl from Grossan, 
12 AGNl from Boyle et al. (1998) and 25 AGNl from Akiyama et al. 
(2000). The excess of faint AGNl at high redshifts is even more evident 
than in the previous figure. The dotted line is the PLE fit (model 0, 
Table 1) by Boyle et al. (1998). The dashed hue is our best PLE fit 
(model 1, Table 1). b) The luminosity function of AGNl fitted with 
the LDDE model (model 4 in Table 1, see text). 



parameters free to vary (model 4, see Figure 2b). In this case the 2DKS test 
gives a probability of 0.47. 

Although we find a better fit of the data with the LDDE model, the ex- 
isting sample is not able to distinguish between PLE and LDDE. In Figure 3a 
the differential counts of the 156 AGNl used in the evaluation of the evolu- 
tion of the LF is shown. The continuous line is the prediction of the LDDE/b 
model while the dashed line is the prediction of the PLE model with Zcut- 
(ilm, riA)=(l. 0,0.0) is assumed. The two models differentiates at fluxes fainter 
than /2_ioA;eF~10~^^er3 cm~'^s~^, where the LDDE/b model predicts a higher 
density of AGNl. The new upcoming fainter surveys from Chandra and XMM- 
Newton will easily test which model is correct. 

Our best fit LF LDDE/b model predicts a contribution of AGNl to the 
2-10 kcV X-ray background of /2-10 = 1-04 x 10~^^ erg cm~^ deg~^, while 
the PLE model with Zcut predicts /2-10 = 0.76 x 10~^^ erg cm~^ s~^ deg~^. 
These values correspond to a fraction of 53% and 39% respectively of the 2-10 
keV X-ray background (we used /2-10 = 1-95 x 10~^^ erg cm~^ s"^ deg~^ from 
Chen, Fabian and Gendreau 1997). 

If an {Qt^,Q,\)=(0.3,0.7) cosmology is assumed the PLE models obtain an 
ever better representation of the data in comparison with what found in the 
(Qm5f^A) = (l-0,0.0) Universe (see Table 1, models 5 and 6). The 2DKS proba- 
bility is 0.73 and 0.48, with and without the introduction of the Zcut parameter 
respectively. In this case even the simple PLE model obtains a quite good fit of 
the data (Figure 3b), and the introduction of the Zcut parameter is necessary to 
stop the evolution only at redshifts larger than 2.3. However our data contain 
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Figure 3. a) The differential counts of AGNl. The continuous hne 
is the expected counts for the LDDE/b model. The dashed line is the 
expected counts for the PLE model with Zcut- i^m, ^^a) = (1-0,0.0) is as- 
sumed, b) The luminosity function in the {O,^,^l\)={0.3,0.7) Universe 
fitted by our PLE model (model 5 in Table 1). 



not enough AGNl at redshift larger than 2 in order to obtain an accurate mea- 
sure of the Zcut parameter (see Figure la), and the errors in this parameter are 
quite large. 



3. ELAIS: The evolution of AGNl in the MIR (15/xm) 
3.1. The ELAIS survey 

The European Large Area ISO Survey (ELAIS) is the largest single open time 
project conduced by ISO, mapping an area of 12 deg^ at 15;um with ISOCAM 
and at 90/xm with ISOPHOT. Four main fields were chosen (Nl, N2, N3 in the 
north hemisphere and SI in the south) at high Echptic latitudes (|/3| > 40°). 



Table 1. The 2-10 keV AGN LF parameters 



Model 


71 


72 


LogL* 


k or pi 


^cut 


A" 








P2DKS 


0) PLE from Boyle* 


1.73 


2.96 


44.16 


2.00 




8.2 


X 


10" 


-r 


0.04" 


1) PLE* 


1.83 


3.00 


44.17 


2.12 




9.8 


X 


10" 


-7 


0.22 


2) PLE with Zcut* 


1.87 


3.03 


44.17 


2.52 


1.39 


8.9 


X 


10" 


-7 


0.31 


3) LDDE/a * 


0.62= 


2.25= 


44.13 


5.40= 


1.55= 


1.4 


X 


10" 


-6 c,d 


0.11 


4) LDDE/b* 

5) PLE* 


0.68 


2.05 


44.03 


4.66 


1.55= 


2.0 


X 


10" 


-6 d 


0.47 


1.93 


2.96 


44.23 


2.19 




8.6 


X 


10" 


-7 


0.48 


6) PLE with Zcut" 


1.92 


2.96 


44.22 


2.23 


2.31 


8.5 


X 


10" 


-7 


0.73 



68% roufidoucn errors i;;;^^ i;;;}:; s'a. 

(*) (Un„UA)=(1.0,0.0j; (o) (i2n„S2A)=(0.3,Q.7j 

(a) Mpc"^(10''^erg s"^)~^ (b) probability; (c) fixed; (d) Mpc-^(erg s"^)"^ 
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Figure 4. a) Integral Counts at 15jUm from the Final Analysis of the 
ELAIS-Sl Region (light grey shaded area). Also plotted the faintest 
IRAS surveys and the deep ISOCAM surveys (dark grey shaded areas), 
b) The luminosity/redshift distribution of our 21 ELAIS-Sl and 41 
RMS AGNls that entered into the computation. 



An initial catalog for SI (J2000,q: : 00^34"'W ,6 : -43°34'44"), covering 
and area of 3.96 deg^, was produced using the Imperial College data reduction 
technique ("Preliminary Analysis", Serjeant et al., 2000). Optical identifications 
were possible thanks to an extensive R-band CCD survey, performed with the 
ESO/Danish 1.5m telescope. The spectroscopic follow-up program, carried out 
at the AAT at the AAO and the 3.6m/NTT at ESO/La Silla, of 114 sources in 
SI fainter than R~17.0 provided the sample presented here. 

The ELAIS-Sl field has also been completely covered in the radio at 1.4 GHz 
down to 0.3 mJy (Gruppioni et al., 1999), and 50% covered in the X-rays with 
BeppoSAX (Alexander et al., 2001). Now the Final Analysis has been completed 
by the ELAIS team in Bologna. The resulting complete catalog includes more 
that 450 sources, with fiuxes at 15/xm down to 0.5 mJy and selected with S/N>5 
(Lari et al., 2001). 

For many years there has been a gap in fiux between the brighter samples 
(>300 mJy), coming from the IRAS surveys and covering large sky areas, and 
the deep/pencil-bin surveys carried out by ISO at much fainter fluxes (<1 mJy). 
In the integral counts derived from the final analysis (Figure 4a, Gruppioni et 
al., in prep) we see how the ELAIS survey fills the whole flux range between 
these two regimes. 

The nature of the objects identified by the spectroscopic follow-up is: 

— High star-forming galaxies (45%); 

— AGNs (type 1 & 2) represent 30%; 

— 15% of galaxies dominated by absorption lines; 

— 4% of late-type stars; 

— 6% of unclear classification (AGN2, Starburst, LINER). 
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Figure 5. a) Distribution of AGNl (plus symbols) and AGN2 

(crossed squares) in the redshift-Rmag space. While the R-magnitude 
distribution is constant with a large spread for AGNl, AGN2 show a 
trend with redshift, becoming fainter than R=20.0 at z ~0.5. b) Dis- 
tribution in the redshift space of AGNl and AGN2 (shaded area). If we 
consider the redshift bin 2;=[0.0,0.5] the ratio of AGN2/AGN1 becomes 
ll/3~4. 



Table 2. The percentages of all objects with emission lines 







AGNl 


AGN2 


ELG 


LINERS 


ELAIS-Sl 


% 


20 


10 


50 


2 


15/im 


<z> 


1.4 


0.3 


0.2 


0.3 


IRAS 


% 


6 


17 


54 


23 


12/im 


<z> 


0.04 


0.015 


0.014 


0.007 


Note - IRAS 


subsample from Alexander & Aussel (2000) 





The observed ratio of the number of type 2 with respect to type 1+2 AGNs is 
1/3. But this result is mainly an artifact due to the different selection functions 
of the sample for AGNl and AGN2. The optical selection introduced in the 
spectroscopic follow-up (17.0<i?<20.0) seems to be at the origin of the lack of 
type 2 AGNs beyond z ^0.5. As a matter of fact, in Figure 5a a trend is observed 
for AGN2 with optical magnitudes being fainter at larger redshifts. In Figure 
5b the redshift distributions of AGNl and AGN2 are shown. In the redshift bin 
where type 2 AGNs are observed (2;=0.0-0.5), the ratio AGN2/AGN1 becomes 
~4, similar to the predictions of the standard unification model. The follow- 
up campaign to be carried out at September /November 2001 with the ESQ 
telescopes will verify if this trend is seen at fainter and brighter optical fluxes. 

3.2. The Evolution of AGNl 

We have assumed Ho=75 km/s/Mpc and a (r2m5^A)=(l-0,0.0) cosmology. Our 
ELAIS preliminary sample of 21 AGNl is statistically significant enough to 
compute a first estimate on the evolution of these objects in the Mid-IR. The 
mean SED from Elvis et al. (1994) for radio quiet QSOs was adopted as a good 
representation of our sources. ivT-corrections were computed for each of the two 
different filters: IS0CAM-LW3 at 15//m and IRAS 12/xm. 

A subsample of AGNl was extracted from the catalog of Rush, Malkan 
&; Spinoglio (1993), as representative in the local universe of this type of ob- 
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jects. The catalog consists of a sample of galaxies selected at 12/xm from the 
IRAS Point Source Catalog PSCv2, and is complete down to 0.3 Jy. With the 
computed X-correction, 15;um uL^ luminosities (L15) were derived. Figure 4b 
represents, in the luminosity-redshift space, all AGNl coming from ELAIS and 
RMS that have been used in this analysis. 

Similarly to what found in the optical and in the X-rays we adopted a 
smooth double power-law for the space density distribution of QSOs in the local 
universe {z=0): 

— = U (LiR/L*j + (LiR/L*j^ , 

dLogLiR L J 

the standard luminosity evolution (PLE) has been adopted: Li^(z) = 
Li5(0)(l+2:)^. A parametric, unbinned maximum likelihood method was used to 
fit the evolution and luminosity function parameters simultaneously (Marshall 
et al., 1983) at Ibfim. Since ELAIS identifications were not only flux limited at 
15//m, but also in their R-band magnitude (17.0 < R < 20.0), a factor was 
introduced in the function 'S' to be minimized 



S = -2 ^ ln[$(zi, Li)] + JJ $(z, L)n(z, L) © "^dzdL, 

to correct for incompleteness, and only applied to the ELAIS sample. This factor 
represents the probability that a source with a given luminosity at 15/im (L15) 
has a R-magnitude (Lr) between the limits of the sample (17.0<R<20.0), 

0(z,L) (17.2 < R < 20.0 | L15) . 

and was derived taking into account the la internal spread in the SED. 

Prom the total available list of AGNl of the sample, the number of objects 
that entered the computation was defined as: 

- RMS: sources with Fi2;um > 300 mJy, 41 sources 

- ELAIS: sources with Fisyj^m > 1 mJy, and 17.0 < R < 20.0, 21 sources 

The resulting estimate of the LP of AGNl at 15/um with its parameters 
is shown in Figure 6. The probability that our data have been drawn from 
the fitted PLE model is 0.28, as given by the 2DKS test. The PLE model 
adopted here could not be sufficient to represent the space density of our sources, 
some density evolution may be required. Currently we are studying different 
parameterizations taking into accoTint some degree of density evolution. 

For the derived PLE model the contribution of AGNl to the CIRB at 
Ibfim is vli, = 5.2 x 10~^^ Wm^^sr^^ which corresponds ~ 2% of the lower 
limit to the CIRB calculated by Altieri et al. 1999 (i//^ = 3.3 x lO'^ W m'^ 
sr~^) in this band (see also Hauser Sz Dwek 2001; Hauser 2001). At maximum, 
the total contribution of AGNs to the background at 15/im can be as high as 
vli, = 2.6 X 10~^ Wm~^sr~^ (~10% of the background measured by Altieri et 
al. under the extreme assumptions: 1) that AGN2 are as much bright as AGNl 
at 15 /um, 2) that the ratio of type 2 to type 1 AGN is 4 at all redshifts, and 2) 
that AGN2s evolved with the same LF of AGNl. 
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Figure 6. PLE fit to our 63 total sources (RMS + ELAIS-Sl). The 
densities have been corrected for evolution within the redshift bins. 



3.3. Conclusions 

Thanks to the spectroscopic campaigns carried out on the HELLAS and ELAIS 

catalogues, we have been able to build up statistical significant samples of AGNl 

and estimate the evolution of their LF in the hard X-rays and IR. 

In the hard X-rays we have been able to double the number of hard X- 

ray AGNl available for statistical analysis at fluxes in the range f2-iokeV ~ 
10-^2 _ 10-13.5 

erg cm s . In total we have used 74 AGNl at these fluxes, 
which combined with the local sample of Grossan (1992) have allowed to show 
directly the shape of the LF of AGNl as function of redshift. In the IR, our 
sample of 21 AGNl has allowed the first estimate of the evolution of the LF at 
15//m. 

Both LFs are fairly well represented by a double-power-law-function with 
a significant cosmological evolution according to a PLE model with L(z)cx(l -|- 
z)*^. However, there are some evidences of density evolution (LDDE) at faint 
luminosities. 
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